Combined pesticides in field doses weaken honey bee (Apis cerana F.) flight ability and analyses of transcriptomics and metabolomics.

Imidacloprid, chlorpyrifos, and glyphosate rank among the most extensively employed pesticides worldwide. The effects of these pesticides and their combined on the flight capability of Apis cerana, and the potential underlying mechanisms remain uncertain. To investigate these effects, we carried out flight mill, transcriptome, and metabolome experiments. Our findings reveal that individual acute oral treatments with pesticides, specifically 20 µL of 10 ng/g imidacloprid (0.2 ng per bee), 30 ng/g chlorpyrifos (0.6 ng per bee), and 60 ng/g glyphosate (1.2 ng per bee), did not impact the flight capability of the bees. However, when bees were exposed to a combination of two or three pesticides, a notable reduction in flight duration and distance was observed. In the transcriptomic and metabolomic analyses, we identified 307 transcripts and 17 metabolites that exhibited differential expression following exposure to combined pesticides, primarily associated with metabolic pathways involved in energy regulation. Our results illuminate the intricate effects and potential hazards posed by combined pesticide exposures on bee behavior. These findings offer valuable insights into the synergistic potential of pesticide combinations and their capacity to impair bee behavior. Understanding these complex interactions is essential for comprehending the broader consequences of pesticide formulations on honey bee populations.

Phylogenomic analysis of Stylops reveals the evolutionary history of a Holarctic Strepsiptera radiation parasitizing wild bees.

Holarctic Stylops is the largest genus of the enigmatic insect order Strepsiptera, twisted winged parasites. Members of Stylops are obligate endoparasites of Andrena mining bees and exhibit extreme sexual dimorphism typical of Strepsiptera. So far, molecular studies on Stylops have focused on questions on species delimitation. Here, we utilize the power of whole genome sequencing to infer the phylogeny of this morphologically challenging genus from thousands of loci. We use a species tree method, concatenated maximum likelihood analysis and Bayesian analysis with a relaxed clock model to reconstruct the phylogeny of 46 Stylops species, estimate divergence times, evaluate topological consistency across methods and infer the root position. Furthermore, the biogeographical history and coevolutionary patterns with host species are assessed. All methods recovered a well resolved topology with close to all nodes maximally supported and only a handful of minor topological variations. Based on the result, we find that included species can be divided into 12 species groups, seven of them including only Palaearctic species, three Nearctic and two were geographically mixed. We find a strongly supported root position between a clade formed by the spreta, thwaitesi and gwynanae species groups and the remaining species and that the sister group of Stylops is Eurystylops or Eurystylops + Kinzelbachus. Our results indicate that Stylops originated in the Western Palaearctic or Western Palaearctic and Nearctic in the early Neogene or late Paleogene, with four independent dispersal events to the Nearctic. Cophylogenetic analyses indicate that the diversification of Stylops has been shaped by both significant coevolution with the mining bee hosts and host-shifting. The well resolved and strongly supported phylogeny will provide a valuable phylogenetic basis for further studies into the fascinating world of Strepsipterans.

Identification of peptides from honeybee gut symbionts as potential antimicrobial agents against Melissococcus plutonius.

Eusocial pollinators are crucial elements in global agriculture. The honeybees and bumblebees are associated with a simple yet host-restricted gut community, which protect the hosts against pathogen infections. Recent genome mining has led to the discovery of biosynthesis pathways of bioactive natural products mediating microbe-microbe interactions from the gut microbiota. Here, we investigate the diversity of biosynthetic gene clusters in the bee gut microbiota by analyzing 477 genomes from cultivated bacteria and metagenome-assembled genomes. We identify 744 biosynthetic gene clusters (BGCs) covering multiple chemical classes. While gene clusters for the post-translationally modified peptides are widely distributed in the bee guts, the distribution of the BGC classes varies significantly in different bee species among geographic locations, which is attributed to the strain-level variation of bee gut members in the chemical repertoire. Interestingly, we find that Gilliamella strains possessing a thiopeptide-like BGC show potent activity against the pathogenic Melissococcus plutonius. The spectrometry-guided genome mining reveals a RiPP-encoding BGC from Gilliamella with a 10 amino acid-long core peptide exhibiting antibacterial potentials. This study illustrates the widespread small-molecule-encoding BGCs in the bee gut symbionts and provides insights into the bacteria-derived natural products as potential antimicrobial agents against pathogenic infections.

Genetic variability of the honey bee mite, Varroa destructor, from a humid continental climatic region of Canada, and temperate and tropical climatic regions of Mexico.

Varroa destructor is a damaging mite of Western honey bees (Apis mellifera). Genetic variability of the mite in different regions of the world could be related to the movement of infested bees or other factors, such as climate. In this study, V. destructor samples were collected from tropical and temperate climate regions of Mexico, and a humid continental climate region of Canada. COX-1 AFLPs showed that all the mites were the Korean haplotype. Four microsatellites revealed nine haplogroups from the continental climate region of Canada, compared to three haplogroups from the tropical and temperate climate regions of Mexico. CytII-ATP sequences showed seven haplogroups from the humid continental climate region vs. two haplogroups from the temperate region and one haplogroup from the tropical region. CytB sequences revealed seven haplogroups from Canada vs. three from Mexico. A comparison of the cytB sequences of the samples from Canada and Mexico to those from a worldwide collection showed that one sequence, designated the cytB1 type, predominated, comprising 57% of the 86 sequences; it clustered with similar sequences that comprised 80% of the sequences, designated family A. CytB1 was predominant in Mexico, but not in Canada. The other 20% of sequences were in families B and C, and all those samples originated from East and Southeast Asia. The microsatellite, cytII-ATP, and cytB markers, all showed higher variability in mites collected in Canada than in Mexico, which could be related to the cooler climate or an earlier invasion and/or multiple mite invasions in Canada.

Ame-miR-980-3p participates in autophagy-mediated midgut remodelling in Apis mellifera via targeting Atg2B.

Autophagy is a process that serves to degrade damaged proteins and organelles, thereby promoting cell homeostasis, differentiation, development and survival. Many miRNAs have been found to have regulatory roles in autophagy. In insects, it has been shown that autophagy is involved in hormone-regulated programmed cell death during metamorphic midgut remodelling. However, whether this is also true during the remodelling of the honey bee midgut is unclear. In the present study, we explored the relationship between autophagy and midgut remodelling and sought to identify miRNAs involved in this physiological process. We found that autophagy occurred during midgut remodelling and that the inhibition of autophagy resulted in midgut dysplasia in prepupae. Differentially expressed miRNAs enriched in the autophagy signalling pathway during midgut remodelling were identified by small RNA-seq. Ame-miR-980-3p, which targets the autophagy-related gene Atg2B, was screened out. Furthermore, abnormal expression of ame-miR-980-3p in the pupal stage led to the thinning of the midgut wall of newly emerged bees (NE). When ame-miR-980-3p expression was inhibited, the intestinal villi of NE bees became significantly shorter and sparse, and the lipid signal in the peritrophic matrix of Pb almost disappeared, indicating that the adult midgut was underdeveloped and the lipid absorption ability was weakened. Taken together, ame-miR-980-3p targeted Atg2B to participate in the regulation of midgut autophagy in the pupae, and the abnormal expression of ame-miR-980-3p would interfere with cell proliferation and death in the process of midgut remodelling, hinder the formation of adult midgut and eventually lead to adult midgut dysplasia and affect the lipid absorption function of the midgut in Apis mellifera.

Apis mellifera Solinvivirus-1, a Novel Honey Bee Virus That Remained Undetected for over a Decade, Is Widespread in the USA.

A metagenomic analysis of the virome of honey bees (Apis mellifera) from an apiary with high rates of unexplained colony losses identified a novel RNA virus. The virus, which was named Apis mellifera solinvivirus 1 (AmSV1), contains a 10.6 kb positive-strand genomic RNA with a single ORF coding for a polyprotein with the protease, helicase, and RNA-dependent RNA polymerase domains, as well as a single jelly-roll structural protein domain, showing highest similarity with viruses in the family Solinviviridae. The injection of honey bee pupae with AmSV1 preparation showed an increase in virus titer and the accumulation of the negative-strand of AmSV1 RNA 3 days after injection, indicating the replication of AmSV1. In the infected worker bees, AmSV1 was present in heads, thoraxes, and abdomens, indicating that this virus causes systemic infection. An analysis of the geographic and historic distribution of AmSV1, using over 900 apiary samples collected across the United States, showed AmSV1 presence since at least 2010. In the year 2021, AmSV1 was detected in 10.45% of apiaries (95%CI: 8.41-12.79%), mostly sampled in June and July in Northwestern and Northeastern United States. The diagnostic methods and information on the AmSV1 distribution will be used to investigate the connection of AmSV1 to honey bee colony losses.

AccsHSP21.7 enhances the antioxidant capacity of Apis cerana cerana.

BACKGROUND: The widespread use of glyphosate has many adverse effects on Apis cerana cerana. Due to the incomplete understanding of the molecular mechanisms of glyphosate toxicity, there are no available methods for mitigating the threat of glyphosate to Apis cerana cerana. Small heat shock proteins (sHSPs) play an important role in resisting oxidative stress, but their mechanism of action in Apis cerana cerana remains unclear. RESULTS: In this experiment, we cloned and identified AccsHSP21.7. Studies have shown that AccsHSP21.7 contains binding motifs for various transcription factors related to oxidative stress. Abiotic stresses induced the expression of AccsHSP21.7. Bacteriostatic testing of a recombinant AccsHSP21.7 protein proved that Escherichia coli overexpressing AccsHSP21.7 showed increased resistance to oxidative stress. Knocking down the AccsHSP21.7 gene caused significant damage to midgut cells, which seriously disrupted the antioxidant system in Apis cerana cerana and greatly increased mortality under glyphosate stress. CONCLUSION: This study investigated the relationship between antioxidant regulation and the AccsHSP21.7 gene at the molecular level, and the results have guiding significance for the improvement of stress resistance in Apis cerana cerana. © 2023 Society of Chemical Industry.

Different effects of pesticides on transcripts of the endocrine regulation and energy metabolism in honeybee foragers from different colonies.

Honeybees are important pollinators of many crops and contribute to biological biodiversity. For years, a decline in bee populations has been observed in certain areas. This decline in honeybees is accompanied by a decrease in pollinator services. One factor contributing to the decline of bee colonies is the exposure to pesticides. Pesticide exposure of bees, among other effects, can negatively affect orientation, memory, immune system function and gene expression. Among the altered expressed genes are transcripts of endocrine regulation and oxidative phosphorylation. Endocrine regulation plays an important role in the development of nurse bees into foragers and oxidative phosphorylation is involved in energy metabolism. Most of these transcriptional changes were investigated using mixed aged honeybees derived from the same colony. Experiments using nurse bees or foragers of the same age but from different colonies are rare. In the present study, effects of the two pesticides chlorpyrifos and pyraclostrobin on the expression of transcripts linked to endocrine regulation and oxidative phosphorylation in foragers of the same age from three different colonies are investigated to fill this gap. These two pesticides were selected because negative effects at sublethal concentrations on bees are known and because they are found in pollen and nectar of crops and wild plants. For this purpose, 20-22 days old foragers of three different colonies were exposed to different sublethal concentrations of the selected fungicides for 24 h, followed by analysis of the expression of buffy, vitellogenin, hbg-3, ilp-1, mrjp1, 2 and 3, cox5a, cox5b and cox17. Some significant changes in gene expression of both endocrine regulation transcripts and oxidative phosphorylation were shown. Furthermore, it became clear that forager bees from different colonies react differently. This is especially important in relation to the risk analysis of pesticides. In addition, it could be shown that the expression of hbg-3 in the brain of bees is a robust marker to distinguish nurse bees from foragers at the molecular biological level. In summary, this study clearly shows that pesticides, which are often detected in pollen and nectar, display negative effects at sublethal concentrations on bees and that it is important to use bees from different colonies for risk assessment of pesticides.

Single gene mutation in a plant MYB transcription factor causes a major shift in pollinator preference.

Understanding the molecular basis of reproductive isolation and speciation is a key goal of evolutionary genetics. In the South American genus Petunia, the R2R3-MYB transcription factor MYB-FL regulates the biosynthesis of UV-absorbing flavonol pigments, a major determinant of pollinator preference. MYB-FL is highly expressed in the hawkmoth-pollinated P. axillaris, but independent losses of its activity in sister taxa P. secreta and P. exserta led to UV-reflective flowers and associated pollinator shifts in each lineage (bees and hummingbirds, respectively). We created a myb-fl CRISPR mutant in P. axillaris and studied the effect of this single gene on innate pollinator preference. The mutation strongly reduced the expression of the two key flavonol-related biosynthetic genes but only affected the expression of few other genes. The mutant flowers were UV reflective as expected but additionally contained low levels of visible anthocyanin pigments. Hawkmoths strongly preferred the wild-type P. axillaris over the myb-fl mutant, whereas both social and solitary bee preference depended on the level of visible color of the mutants. MYB-FL, with its specific expression pattern, small number of target genes, and key position at the nexus of flavonol and anthocyanin biosynthetic pathways, provides a striking example of evolution by single mutations of large phenotypic effect.

Chromatin accessibility-based characterisation of brain gene regulatory networks in three distinct honey bee polyphenisms.

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

Dual function of a bumblebee (Bombus ignitus) serine protease inhibitor that acts as a microbicidal peptide and anti-fibrinolytic venom toxin.

In bee venoms, low-molecular-weight peptides, including serine protease inhibitors (SPIs), exhibit multifunctional activities. Although SPIs in bee venoms are relatively well known, those that function in both the body and secreted venom of bees are not well-characterized. In this study, we identified a bumblebee (Bombus ignitus) SPI (BiSPI) that displays microbicidal and anti-fibrinolytic activities. BiSPI was found to consist of a trypsin inhibitor-like domain containing a P1 site and ten cysteine residues. We observed that the BiSPI gene was ubiquitously transcribed in the body, including the venom glands. In correlation, the BiSPI protein was detected both in the body and secreted venom by using an antibody against a recombinant BiSPI peptide produced in baculovirus-infected insect cells. Recombinant BiSPI exhibited inhibitory activity against trypsin but not chymotrypsin and inhibited microbial serine proteases and plasmin but not elastase or thrombin. Moreover, recombinant BiSPI recognized carbohydrates and bound to fungi and gram-negative and gram-positive bacteria. Consistent with these properties, recombinant BiSPI exhibited microbicidal activities against bacteria and fungi through induction of structural damage by binding to the microbial surfaces. Additionally, recombinant BiSPI inhibited the plasmin-mediated degradation of human fibrin and was thus concluded to exhibit anti-fibrinolytic activity. Moreover, the peptide showed no effect on hemolysis. These findings demonstrate the dual function of BiSPI, which acts as a microbicidal peptide and anti-fibrinolytic venom toxin.

DWV 3C Protease Uncovers the Diverse Catalytic Triad in Insect RNA Viruses.

Deformed wing virus (DWV) is the most prevalent Iflavirus that is infecting honey bees worldwide. However, the mechanisms of its infection and replication in host cells are poorly understood. In this study, we analyzed the structure and function of DWV 3C protease (3Cpro), which is necessary for the cleavage of the polyprotein to synthesize mature viral proteins. Thus, it is one of the nonstructural viral proteins essential for the replication. We found that the 3Cpros of DWV and picornaviruses share common enzymatic properties, including sensitivity to the same inhibitors, such as rupintrivir. The predicted structure of DWV 3Cpro by AlphaFold2, the predicted rupintrivir binding domain, and the protease activities of mutant proteins revealed that it has a Cys-His-Asn catalytic triad. Moreover, 3Cpros of other Iflaviruses and Dicistrovirus appear to contain Asn, Ser, Asp, or Glu as the third residue of the catalytic triad, suggesting diversity in insect RNA viruses. Both precursor 3Cpro with RNA-dependent RNA polymerase and mature 3Cpro are present in DWV-infected cells, suggesting that they may have different enzymatic properties and functions. DWV 3Cpro is the first 3Cpro characterized among insect RNA viruses, and our study uncovered both the common and unique characteristics among 3Cpros of Picornavirales. Furthermore, it would be possible to use the specific inhibitors of DWV 3Cpro to control DWV infection in honey bees in future. IMPORTANCE The number of managed honey bee (Apis mellifera) colonies has considerably declined in many developed countries in the recent years. Deformed wing virus (DWV) vectored by the mites is the major threat to honey bee colonies and health. To give insight into the mechanism of DWV replication in the host cells, we studied the structure-function relationship of 3C protease (3Cpro), which is necessary to cleave a viral polyprotein at the specific sites to produce the mature proteins. We found that the overall structure, some inhibitors, and processing of 3Cpro are shared between Picornavirales; however, there is diversity in the catalytic triad. DWV 3Cpro is the first viral protease characterized among insect RNA viruses and reveals the evolutionary history of 3Cpro among Picornavirales. Furthermore, DWV 3Cpro inhibitors identified in our study could also be applied to control DWV in honey bees in future.

Molecular phylogeny, historical biogeography and revised classification of andrenine bees (Hymenoptera: Andrenidae).

The mining bee subfamily Andreninae (Hymenoptera: Andrenidae) is a widely distributed and diverse group of ground-nesting solitary bees, including numerous species known to be important pollinators. Most of the species diversity of Andreninae is concentrated in the mainly Holarctic genus Andrena, comprising ca. 1550 described species. The subfamily and especially the genus have remained relatively neglected by recent molecular phylogenetic studies, with current classifications relying largely on morphological characters. We sampled ultraconserved element (UCE) sequences from 235 taxa, including all andrenine genera and 98 out of 104 currently recognized Andrena subgenera. Using 419,858 aligned nucleotide sites from 1009 UCE loci, we present a comprehensive molecular phylogenetic analysis of the subfamily. Our analysis supports the recognition of seven distinct genera in the Andreninae: Alocandrena, Ancylandrena, Andrena, Cubiandrena, Euherbstia, Megandrena, and Orphana. Within the genus Andrena, present-day subgeneric concepts revealed high degrees of paraphyly and polyphyly, due to strong homoplasy of morphological characters, necessitating a thorough, extensive revision of the higher classification of the genus. Based on our findings, we place the subgenus Calcarandrena in synonymy with Andrena (Lepidandrena); Hyperandrena, Nemandrena, Scoliandrena, Tylandrena and Zonandrena with A. (Melandrena); Distandrena, Fumandrena and Proxiandrena with A. (Micrandrena); Carandrena with A. (Notandrena); Agandrena with A. (Plastandrena); Geandrena and Xanthandrena with A. (Ptilandrena); Xiphandrena with A. (Scrapteropsis); and Platygalandrena and Poliandrena with A. (Ulandrena) (new synonymies). We additionally reestablish the groups known as Opandrena and Truncandrena as valid subgenera of Andrena. Our results also show that the MRCA of Andrena + Cubiandrena dispersed from the New World to the Palaearctic probably during the Eocene-early Oligocene, followed by 10-14 Neogene dispersal events from the Palaearctic to the Nearctic and 1-6 Neogene dispersals back into the Palaearctic, all within the genus Andrena.

Phylogenetic relationships and the evolution of host preferences in the largest clade of brood parasitic bees (Apidae: Nomadinae).

Brood parasites (also known as cleptoparasites) represent a substantial fraction of global bee diversity. Rather than constructing their own nests, these species instead invade those of host bees to lay their eggs. Larvae then hatch and consume the food provisions intended for the host's offspring. While this life history strategy has evolved numerous times across the phylogeny of bees, the oldest and most speciose parasitic clade is the subfamily Nomadinae (Apidae). However, the phylogenetic relationships among brood parasitic apids both within and outside the Nomadinae have not been fully resolved. Here, we present new findings on the phylogeny of this diverse group of brood parasites based on ultraconserved element (UCE) sequence data and extensive taxon sampling with 114 nomadine species representing all tribes. We suggest a broader definition of the subfamily Nomadinae to describe a clade that includes almost all parasitic members of the family Apidae. The tribe Melectini forms the sister group to all other Nomadinae, while the remainder of the subfamily is composed of two sister clades: a "nomadine line" representing the former Nomadinae sensu stricto, and an "ericrocidine line" that unites several mostly Neotropical lineages. We find the tribe Osirini Handlirsch to be polyphyletic, and divide it into three lineages, including the newly described Parepeolini trib. nov. In addition to our taxonomic findings, we use our phylogeny to explore the evolution of different modes of parasitism, detecting two independent transitions from closed-cell to open-cell parasitism. Finally, we examine how nomadine host-parasite associations have evolved over time. In support of Emery's rule, which suggests close relationships between hosts and parasites, we confirm that the earliest nomadines were parasites of their close free-living relatives within the family Apidae, but that over time their host range broadened to include more distantly related hosts spanning the diversity of bees. This expanded breadth of host taxa may also be associated with the transition to open-cell parasitism.

Telomerase-Positive Somatic Tissues of Honeybee Queens (Apis mellifera) Display No DNA Replication.

Telomere biology is closely linked to the process of aging. The restoration of telomere length by maintaining telome-rase activity in certain cell types of human adults allows for the proliferative capacity of the cells and preserves the regeneration potential of the tissue. The absence of telome-rase, that leads to telomere attrition and irreversible cell cycle arrest in most somatic cells, acts as a protective mechanism against uncontrolled cancer growth. Nevertheless, there have been numerous studies indicating noncanonical functions of telomerase besides those involved in telomere lengthening. Eusocial insects serve as a great system for aging research. This is because eusocial reproductives, such as queens and kings, have a significantly extended lifespan compared to nonreproductive individuals of the same species. We report that the somatic tissues of honeybee queens (Apis mellifera) are associated with upregulated telomerase activity; however, this upregulation does not fully correlate with the rate of DNA replication in the tissues. This indicates a noncanonical role of telomerase in the somatic tissues of honeybee queens.

AccPDIA6 from Apis cerana cerana plays important roles in antioxidation.

PDIA6 is a member of the protein disulfide isomerase (PDI) family, shows disulfide isomerase activity and oxidoreductase activity, and can act as a molecular chaperone. Its biological functions include modulating apoptosis, regulating the proliferation and invasion of cancer cells, supporting thrombosis and modulating insulin secretion. However, the roles of PDIA6 in Apis cerana cerana are poorly understood. Herein, we obtained the PDIA6 gene from A. cerana cerana (AccPDIA6). We investigated the expression patterns of AccPDIA6 in response to oxidative stress induced by H2O2, UV, HgCl2, extreme temperatures (4 °C, 42 °C) and pesticides (thiamethoxam and hexythiazox) and found that AccPDIA6 was upregulated by these treatments. Western blot analysis indicated that AccPDIA6 was also upregulated by oxidative stress at the protein level. In addition, a survival test demonstrated that the survival rate of E. coli cells expressing AccPDIA6 increased under oxidative stress, suggesting a possible antioxidant function of AccPDIA6. In addition, we tested the transcripts of other antioxidant genes and found that some of them were downregulated in AccPDIA6 knockdown samples. It was also discovered that the antioxidant enzymatic activity of superoxide dismutase (SOD) decreased in AccPDIA6-silenced bees. Moreover, the survival rate of AccPDIA6 knockdown bees decreased under oxidative stress, implying that AccPDIA6 may function in the oxidative stress response by enhancing the viability of honeybees. Taken together, these results indicated that AccPDIA6 may play an essential role in counteracting oxidative stress.

Pollinator assemblage and pollen load differences on sympatric diploid and tetraploid cytotypes of the desert-dominant Larrea tridentata.

PREMISE: Whole-genome duplication (polyploidy) is an important force shaping flowering-plant evolution. Ploidy-specific plant-pollinator interactions represent important community-level biotic interactions that can lead to nonrandom mating and the persistence of mixed-ploidy populations. METHODS: At a naturally occurring diploid-tetraploid contact zone of the autopolyploid desert shrub Larrea tridentata, we combined flower phenology analyses, collections of bees on plants of known cytotype, and flow cytometry analyses of bee-collected pollen loads to investigate whether (1) diploid and tetraploid plants have unique bee pollinator assemblages, (2) bee taxa exhibit ploidy-specific visitation and pollen collection biases, and (3) specialist and generalist bee taxa have ploidy-specific visitation and pollen collection biases. RESULTS: Although bee assemblages overlapped, we found significant differences in bee visitation to co-occurring diploids and tetraploids, with the introduced honeybee (Apis mellifera) and one native species (Andrena species 12) more frequently visiting tetraploids. Consistent with bee assemblage differences, we found that diploid pollen was overrepresented among pollen loads on native bees, while pollen loads on A. mellifera did not deviate from the random expectation. However, mismatches between the ploidy of pollen loads and plants were common, consistent with ongoing intercytotype gene flow. CONCLUSIONS: Our data are consistent with cytotype-specific bee visitation and suggest that pollinator behavior contributes to reduced diploid-tetraploid mating. Differences in bee visitation and pollen movement potentially contribute to an easing of minority cytotype exclusion and the facilitation of cytotype co-occurrence.

RNA m6A Modification Functions in Larval Development and Caste Differentiation in Honeybee (Apis mellifera).

Genetically identical female honeybee larvae with different diets develop into sterile workers or fertile queens. It remains unknown whether the reversible RNA N6-methyladenosine (m6A) mark functionally impact this "caste differentiation." Here, we profile the transcriptome-wide m6A methylome of honeybee queen and worker larvae at three instar stages and discover that m6A methylation dynamics are altered by differential feeding. Multiple methylome comparisons show an obvious increase in m6A marks during larval development and reveal a negative correlation between gene expression and m6A methylation. Notably, we find that worker larvae contain more hypermethylated m6A peaks than do queen larvae, and many caste-differentiation-related transcripts are differentially methylated. Chemical suppression of m6A methylation in worker larvae by 3-deazaadenosine (DAA) reduces overall m6A methylation levels and triggers worker larvae to develop queen caste features. Thus, our study demonstrates that m6A functionally impacts caste differentiation and larval development, yet it does not exclude potential contributions from other factors.

Caste-specific gene expression underlying the differential adult brain development in the honeybee Apis mellifera.

Apis mellifera adult workers feature more developed key brain regions than queens, which allows them to cope with the broad range of duties they need to perform in a colony. However, at the end of larval development, the brain of queens is largely more developed than that of workers. Major morphogenetic changes take place after metamorphosis that shift caste-specific brain development. Here, we tested the hypothesis that this phenomenon is hormonally governed and involves differential gene expression. Our molecular screening approach revealed a set of differentially expressed genes in Pp (first pharate-adult phase) brains between castes mainly coding for tissue remodelling and energy-converting proteins (e.g. hex 70a and ATPsynß). An in-depth qPCR analysis of the transcriptional behaviour during pupal and pharate-adult developmental stage in both castes and in response to artificially augmented hormone titres of 18 genes/variants revealed that: i. subtle differences in hormone titres between castes might be responsible for the differential expression of the EcR and insulin/insulin-like signalling (IIS) pathway genes; ii. the morphogenetic activity of the IIS in brain development must be mediated by ILP-2, iii. which together with the tum, mnb and caspase system, can constitute the molecular effectors of the caste-specific opposing brain developmental trajectories.

Methionine as a methyl donor regulates caste differentiation in the European honey bee (Apis mellifera).

Nutrition contributes to honey bee caste differentiation, but the role of individual nutrients is still unclear. Most essential amino acid contents, except that of methionine (Met), are greater in royal jelly than worker jelly. After ∼3.5 d, the Met content in the latter was slightly greater than in the former. Met is the major raw material used in the synthesis of S-adenosyl-L-methionine, an active methyl donor for DNA methylation, which is an epigenetic driver of caste differentiation. Here, we tested whether Met regulates caste differentiation in honey bees by determining its effects on the caste development of bees receiving four diets: the basic, basic + 0.2% Met, basic + 0.2% Met + 20 mg/kg 5-azacytidine, and basic + 20 mg/kg 5-azacytidine. The presence of Met decreased the adult bee body length and the numbers of ovarioles, indicating that Met may direct the development of female larvae toward worker bees. The upregulated expression of SAMS, Dnmt1, and Dnmt3 caused by Met exposure in 4-d-old larvae indicated that the worker-inductive effects of Met may occur through the promotion of DNA methylation. We investigated the co-effects of Met and glucose on bee development, and found that the effects of an increased glucose level on the number of ovarioles and body length did not strengthen the worker-inductive effects caused by Met. Our results contribute to caste development theory and suggest that Met-as a methyl donor-plays a regulatory, but not decisive, role in caste differentiation.